Single atom adhesion in optimized gold nanojunctions 
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We study the interaction between single apex atoms in a metallic contact, using the break junction 
geometry. By carefully 'training' our samples, we create stable junctions in which no further atomic 
reorganization takes place. This allows us to study the relation between the so-called jump out 
of contact (from contact to tunnelling regime) and jump to contact (from tunnelling to contact 
regime) in detail. Our data can be fully understood within a relatively simple elastic model, where 
the elasticity k of the electrodes is the only free parameter. We find 5 < k < 32 N/m. Furthermore, 
the interaction between the two apex atoms on both electrodes, observed as a change of slope in 
the tunnelling regime, is accounted for by the same model. 



Many macroscopic phenomena find their origin on 
the nanoscale, since they are ultimately due to the 
interaction between single atoms. A good example is 
formed by friction and wear, which have been studied 
for centuries, but still inspire fascinating research. For 
example, several groups have recently explored methods 
to minimize friction in nanoelectromechanical systems 
(NEMS), where no liquid lubricants can be applied [1]. 
In this Letter, we focus on the ultimate miniaturization 
of the problem and investigate adhesion and elasticity 
on the atomic scale. 
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FIG. 1: Points: conductance G vs. distance D for four succes- 
sive Co-loops (Vbias — 50 mV). The jump to contact occurs 
at D=1.5 A] the jump out of contact at D=2.0 A. Black line: 
fit to model in Fig. [2] (k=15.7 N/m). For the other param- 
eters we use literature values: Fq — 1.5nN, d = 2.5 A, and 
Eb — 0.7eV. The work function, cj> — 5eV, was measured 
independently. Inset graph: zoom of G vs. D in contact 
regime (linear scale). Picture: scanning electron micrograph 
of a lithographic break junction. 

The interaction between single atoms can be studied 
by carefully extending a notched metallic wire, while 



monitoring its conductance G. As the wire is thinned 
out, G decreases, until its value is dominated by a 
few atoms forming a constriction 0. When pulling is 
continued, an abrupt rupture of the wire is observed 
(see Fig. [TJ point c to d). Upon closing the contacts, 
a second jump occurs for many metals, including gold 
(see Fig. [TJ point a to b ) 0, H|. These jumps are 
known as the 'jump out of contact' (JOC) and 'jump 
to contact' (JC), respectively. By carefully studying 
these discontinuities, one can in principle obtain detailed 
information on the adhesion forces between two single 
atoms. However, the details of the hysteretic loop in 
Fig. [I] are still not fully understood. In fact, no relation 
between the JC and JOC has been observed so far. The 
reason for this is that the breaking process is generally 
accompanied by plastic deformation, i.e., the atoms first 
reorganize before rupture [5]. Therefore, during closing 
and opening, the electrodes have a different atomic 
configuration. The most intriguing example of plastic 
deformation is the formation of atomic chains prior to 
breaking [6j]. 

Here, we explore junctions in which no plastic deforma- 
tion occurs during breaking and making of the single 
atom contact. This is achieved by properly 'training' 
each device first. Figure [1] displays G during the opening 
and closing of a 'trained' Au wire. Note that the single 
gold atom conductance is characterized by a value close 
to the conductance quantum Go = 2e 2 /h. As can 
be seen from the absence of conductance steps in the 
contact regime, no atomic reorganization takes place. 
Moreover, the curves in Fig. [T]are perfectly reproducible 
for tens of subsequent runs. Being able to exclude 
plasticity, we infer that the two jumps (JC and JOC) 
in Fig. [1] are related to the adhesive forces between 
the single atoms forming the junction. The remarkable 
reproducibility of the 'trained' junctions allows us to test 
a generic potential energy model. In fact, we show that 
the whole making and breaking process can be fitted by 
a single fit parameter: the elasticity of the electrodes. 

For our study, it is crucial to minimize drift and 
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FIG. 2: Total energy as a function of the inter- atomic distance 
x of a gold dimer, for different electrode separations D (see 
inset; Do denotes an offset distance). Two contributions are 
included: one due to the springs (spring constant k) and one 
due to the dimer (described by the 'universal' binding curve) . 
Depending on D, the total energy may have two minima. For 
the atom to jump in and out of contact, a barrier has to 
be overcome. Hence, opening and closing occurs at different 
positions, explaining the hysteresis in Go loops. Note: same 
parameters as in Fig. [T] 



vibrations in the electrodes during the measurements. 
This is achieved by using a mechanically controllable 
break junction (MCBJ) at 4.2 K [7]. The wires are 
patterned with electron beam lithography, after which 
we evaporate 1 nm of Cr and 120 nm Au at 2 * 10 -7 
mbar. Finally, the area below the wire is etched with 
a CF4/O2 plasma, to create a free hanging gold wire 
[inset in Fig. CQ. The MCBJ is cooled down to 4.2 K 
and broken by bending the substrate. Thus, clean and 
stable gold electrodes are obtained in cryogenic vacuum. 
Opening and closing of the junction is done with an 
effective speed of 0.5 A/ s, while the conductance, at a 
50 mV bias, is monitored. Measurements at 1-300 mV 
gave similar results [8f. We emphasize the impressive 
stability of the electrodes, resulting in a drift below 
0.3 pm/h. Our break junctions are calibrated using 
Gundlach oscillations 0, Hol|. We find an attenuation 
factor r = (5.4 ± 0.6) • 10~ 5 [3, El, and a work function 
(j) = beV. 

The details of our 'training' method are as follows. 
When closing the electrodes, we stop immediately as 
soon as the electrodes are in contact, i.e. at G w Go, 
preventing further disorder. Subsequently, we break 
the wire, extend it 1-2 A into the tunnelling regime 
and close it again until the jump to contact. Repeating 
this procedure rearranges and orders the tip atoms. In 
this way, the atoms are able to probe (energetically 
and spatially) different positions, allowing them to 
find the most stable configuration. Remarkably, after 
typically >10 sweeps, JC and JOC occur at two exactly 
reproducible positions. In Fig. [TJ four subsequent 
traces are shown with perfect repeatability. In fact, the 



maximum variation in the closing and opening points 
was less than 5 pm over 50 sweeps. Although these 
loops (which we call 'Go -loops') have already been 
observed by other groups 0, IHj], we are the first to 
optimize the training method to investigate JC and 
JOC in well-defined geometries. We have measured 734 
different Go loops, on 8 different samples, to study the 
variation in the contacts. To obtain a new Go-loop, 
we first rearrange a contact by closing up to > 10Go, 
before training the contact for a different Go-loop. 
For each Go-loop, we automatically record the four 
conductance values G a , G^, G c en Gd (at points a,b,c 
and d in Fig. [TJ respectively). The fact that 85% of the 
conductance values G c is above 0.9Go emphasizes the 
good definition of our junctions. Recently, Untiedt et al 
studied JC (no 'training' method was employed) [l4j. A 
statistical analysis was made of many closing traces and 
correlations were found between the conductance values 
just before and just after JC (G a and For gold, 

they observed maxima in density plots, for G b values 
below Go and around 1.6 Go- Conductances below 1 Go 
were attributed to a dimer configuration, whereas higher 
conductances were related to monomer and double bond 
configurations. We have also observed the peak around 
1.6 Go for untrained junctions. However, upon breaking 
G dropped in steps to lower conductances, making it 
impossible to create stable 'Go-loops' for this configura- 
tion. Furthermore, for our trained contacts, more than 
80 % of our Gb values have a conductance below 1.02 
Go- Therefore, we conclude that training of the contacts 
results predominantly in the dimer configuration, as 
sketched in Fig. [TJ Also from molecular dynamics 
simulations, dimer s are expected to be the most stable 
geometry |14l. Il5j. 

Since for the trained contacts all plastic deformation 
is removed, we can employ an elastic model to describe 
the hysteretic loop in Fig. [TJ The basic configuration is 
shown in the inset of Fig. [2j where a dimer is depicted 
in between two elastic electrodes. To describe the force 
between the two atoms, we use the so called 'universal' 
binding curve, which is given by [HI, [l?], EH 

E{x) = -a(x - x )e-^ x - xo) (1) 

where x is the interatomic distance. The parameters 
a, /3 and xo are related to the equilibrium bond distance 
d = xq + 1//3, the binding energy = —a/ '(3e and 
the slope at the inflection point Fo = a/e 2 , with 
e = 2.718. These values are known from literature, i.e., 
d = 2.5 ± 0.2A pij , F = 1.5 ± 0.3niV for the break force 
[H and E b = 0.7 ± 0.2eV for the binding energy 0. 
For the bonding energy of the dimer to the rest of the 
electrodes (the "banks"), we take the potential energy 
of a spring (ku 2 /2). Hence, we are left with only one 
free parameter, the spring constant k, which is directly 
related to the hysteresis of the Go-loop. In Fig. [2j the 
total energy is plotted versus the interatomic distance 
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for different electrode distances D. Depending on D, 
two minima are present. When starting with a closed 
contact in equilibrium, i.e. D=0 A, the equilibrium 
interatomic distance equals x eq — 2. 5 A (minimum of 
curve 1). As the electrodes are pulled apart, e.g., by 1.8 
A (curve 3), the two atoms are separated by only 0.2 
A. The rest of the displacement is invested in stretching 
the spring. Increasing D further (towards curve 4), the 
first minimum disappears and the system jumps to the 
second minimum (at x eq « 5.3 A); this is JOC. The 
atoms of the dimer get separated by typically 2-3 A. 
Upon closing the junction, a linear relation between x eq 
and D is initially seen. However, at small separations, 
the spring stretches somewhat due to the attractive 
forces of the opposing atoms, i.e. x eq moves faster 
than D [l3|, [2(j • This gives a deviation from exponential 
tunnelling, as observed in Fig. [TJ The effect is maximal 
just before JC, which takes place in between curves 3 
and 2, and covers a distance of « 1A. To apply our 
model to the tunnelling part of the Go-loops, we assume 
that the work function, (/>, does not depend on x eqi so 
that G oc exp(—2x eq ^2m(j)/Ti). Taking G = Go for 
D = 0, we have the tools to fit the data in Fig. [TJ The 
corresponding trace is shown in Fig. [TJ It gives a perfect 
fit, not only to the exact position of JC and JOC, 
but also to the deviation from exponential tunnelling. 
For this Go-loop, the fit parameter k assumes a value 
k = 15.7 N/m. 
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FIG. 3: (a) Conductance G just before JC (G a ) vs. G just 
after JOC (Gd) for 734 different Go-loops, (b) Average G a 
as a function of Gd- Averaging is done within regular bins 
of Gd, as indicated in a). The line is a fit to the model, 
assuming a varying spring constant 5 < k < 32 N/m. The 
other parameters are the same as in Fig. [T]and[2] 

In total, we studied 734 Go-loops, which could all be 
fitted with the model. In Fig. [3j the conductance G a 



(just before JC) is plotted versus Gd (just after JOC), 
for all Go-loops measured. Remarkably, the graph shows 
a relation between these points. This is especially visible 
in Fig. [3)o, where the average G a is displayed versus Gd- 
We compare these data points to our model. By varying 
the elasticity k from 5-32 N/m, with all other parameters 
fixed to literature values, one obtains the line drawn in 
Fig. [3j Clearly, the data points are well described by the 
model, using only k as a variable. Note that the range of 
k- values is in agreement with Ref. [13]. Moreover, > 90 
% of our k values (per electrode) are in the range 7-26 
N/m. This variation is substantially smaller than the 
spread found in Ref [l3[ . We relate the spread in k to the 
following phenomena. First, breaking a gold wire is only 
possible along certain crystal orientations ([111], [100] 
and [110]) [2l|. For each orientation, the apex atom is 
bonded differently to the second layer of the electrode. 
In fact, the apex atom has 3, 4 and 5 nearest neighbors 
for the Au [111], [100] and [110] direction, respectively. 
This is expected to have substantial influence on the 
elasticity of the electrode. However, this picture is not 
yet complete. The work by Olesen implies that k is only 
partly determined by the nearest neighbors of the apex 
atom [22]. In fact, the most significant contribution to k 
arises from elastic displacements in the rest of the metal 
tips. Hence, k is related to the precise structure of the 
electrodes on a larger scale, which varies with each Go 
loop. Every time we close the junction up to 10 Go, we 
most likely introduce defects in the atomic layers further 
into the contact [23], which influence the elasticity of 
the electrodes. A final source of variation is the so-called 
"lateral approach" of the electrodes. For the model in 
Fig. [2j we assumed that the apex atoms are perfectly 
aligned. However, small misalignments (0-1 A) may 
occur in reality. We extended the ID model to a 2D 
model by assuming springs in both x and y directions. 
This yields a variation in k of up to 10 %. We stress 
that the other parameters in our model cannot explain 
the data in Fig. [3j Only the spring constant gives a 
relation along the direction shown. The variation in the 
vertical direction (Fig. [3^), however, can be explained 
by small deviations in E^. As indeed shown by Ref. [l6[, 
Eb is sensitive to the local atomic configuration. We 
find that a spread of 25 % in explains the variation 
in G a . Such a spread is consistent with the results of 
Ref. [13]. We conclude that our relatively simple model 
does not only explain the occurrence of JC and JOC, 
but also their relationship. Finally, we emphasize that 
without 'training', atomic reorganization upon opening 
and closing destroys the interdependence between G a 
and Gd- 

Remarkably, our model fits the data in Figs. [1] and [3] 
very well, despite the fact that we assume a constant 
barrier height <p for all x eq . This contrasts computations 
by Lang which predict a strong decrease of the apparent 
barrier at electrode distances < AA [24|]. Such barrier 
lowering would primarily be due to the image forces 
and local effects related to the electric charges on 
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the electrodes. Olesen et al carefully examined the 
tunnel curves of Ni, Pt and Au by scanning tunnelling 
microscopy (STM) and found no deviations from ex- 
ponential behavior. This was explained by assuming 
that barrier lowering is exactly cancelled by adhesion 
between tip and sample [22|. The fact that we can fit 
our Go-loops using an elastic model only, shows that 
the influence of image forces in break junctions (which 
feature two sharp tips) is relatively small. 




FIG. 4: (a) Conductance just before JOC (G c ) vs. G just 
after JC (Gb) for 734 different Co-loops. The dashed line 
represents G c = Gb- Inset: schematic junction; scale bar: 
barrier length for tunnelling. 

Having described three quarters of the Go-loop (JOC, 
tunnelling and JC), we focus on the contact regime. In 
Fig. [H G c is plotted versus G5. In all cases, we have 



G c < Go, as expected for the conductance of a single 
atom contact. In contrast, the values of Gb go well 
above 1 Go- On average, Gb is 0.025 Go higher than 
G c . This is exactly as expected, if a second conductance 
channel is considered. As calculated by Ref. [15], the 
small distance between the second layers of the two 
electrodes allows for tunnelling with a conductance 
of at most 0.03 Go- When stretching the contact, 
however, the tunnel gap increases and the transmis- 
sion of the second channel tends to zero. This is why 
it is not visible in G c , where the contact is fully extended. 

In summary, we have created highly ordered gold elec- 
trodes, using a "training" method. Upon opening and 
limited closing, our junctions show no plastic deforma- 
tion, allowing us to study the jump out of contact, tun- 
nelling curves and jump to contact in detail. Individual 
breaking and making loops can be perfectly fitted with 
an elastic model, having the spring constant of the elec- 
trodes, k, as the only free parameter. Hence, by suppress- 
ing plastic deformation effects, we are able to measure 
and model adhesion on the single atomic level. 
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